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Abstract 

We have investigated the conditions necessary to achieve stronger Cherenkov-like instability 
of plasma waves leading to emission in the terahertz (THz) regime for semiconductor quantum 
wells (QWs). The surface response function is calculated for a bilayer two-dimensional electron gas 
(2DEG) system in the presence of a periodic spatial modulation of the equilibrium electron density. 
The 2DEG layers are coupled to surface plasmons arising from excitations of free carriers in the 
bulk region between the layers. A current is passed through one of the layers and is characterized 
by a drift velocity V£> for the driven electric charge. By means of a surface response function 
formalism, the plasmon dispersion equation is obtained as a function of frequency cj, the in-plane 
wave vector qii = (q x ,Qy) and reciprocal lattice vector nG where n = 0, ±1, ±2, • • • and G = 2ir/d 
with d denoting the period of the density modulation. The dispersion equation, which yields the 
resonant frequencies, is solved numerically in the complex w-plane for real wave vector qii. It is 
ascertained that the imaginary part of u is enhanced with decreasing d, and with increasing the 
doping density of the free carriers in the bulk medium for fixed period of the spatial modulation. 
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I. INTRODUCTION 



The abundance and availability of information in the literature today is both a challenge 
and a blessing to physicists interested in studying sources of terahertz (THz) radiation. One 
THz (= 10 12 Hz) encompasses frequencies invisible to the naked eye in the electromagnetic 
spectrum, lying between microwave and infrared. Gone are the days when experimentalists 
relied heavily on intuition in their quest for stronger and more uniform sources of radi- 
ation. Instead, todays experimentalists are surrounded by a flood of published material 
yearly. Harvesting this abundance of knowledge and using it to inform and improve this 
area of research presents a challenge for researchers. It is imperative for experimentalists 
to review, critique and transfer scientific findings to effective quality outcomes. One active 
field being pursued is that of THz quantum cascade lasers (QCLs). QCLs were invented 
in 1994 1 , are semiconductor lasers obtained by epitaxially growing a sequence of layers of 
different semiconductors. As a result, electrons in the conduction band are subjected to an 
artificial, periodic, one dimensional potential that varies on the nanometer scale. So far, 
the focus has been on the design, growth, fabrication, and testing of high power THz QCLs 
emitting across a broad frequency range. The work reported up to this point covers sev- 
eral ambitious projects from ultra- long wavelength emission, phase/mode- locking, multiple 
color generation, photonic crystal structures, and improved laser performance with respect 
to both maximum temperature operation and peak output power. 

Far-infrared (FIR) radiation has been used to excite plasmon modes resonantly in a 
quantum-well transistor at frequencies between 0.5 and l.OxlO 12 s" 1 . A split grating gate 
design that allows localized pinch-off of the transistor channel has been found to significantly 
enhance FIR response and to allow electrical tuning of the plasmon resonance.- 1 ^ As a 
matter of fact, the role played by plasma excitations in the THz response of low- dimensional 
microstructured semiconductors has been investigated considerably in recent times.- - — This 
paper discusses a specific set of simulations which bolster earlier work which has already 
been published. We carry out a comprehensive study and report our results for investigating 
plasma instabilities in a pair of Coulomb coupled two-dimensional electron gas (2DEG) 
systems in which inter-layer hopping between layers is not included.— 1 ^ 
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II. GENERAL FORMULATION OF THE PROBLEM 



Our system consists of a pair of parallel two-dimensional electron gas layers embedded 
in the half-space z > 0. The space between the layers is filled by a dielectric material with 
dielectric constant e{ui). In Fig. 1 we show a schematic of the bilayer structure which we 
do numerical calculations for. However, in our formulation of the problem, we consider 
an arbitrary number of parallel 2DEG layers. We assume that our multi-layered structure 
interacts with an external time dependent electrostatic potential e xt(r, t). 

We denote the electric field due to an external electrostatic potential by E cxt = — V0 ex t- 
Let the layered system occupy the region z > 0. When we periodically modulate the 
equilibrium electron density, the surface response function will be modified from its value 
in the absence of the modulation.— 1 ^ We assume that the periodic density modulation is 
along the x-direction and that the modulated sheet density can be described by n 2 n(x) = 

oo 

p n exp(inGx), where the real- value p n is the nth Fourier component of ti2d{x), G = 

71= — OO 

2%/dis the reciprocal lattice vector and d is the period of the modulation. Since V 2 ex t = 
for z > 0, the electrostatic potential in this region and in the vicinity of the surface can 
be written as a superposition of plane waves. This follows from the fact that in the region 
outside the surface, where there is no charge present, the external potential ex t satisfies 
V 2 0cxt = 0. Since the system is translationally invariant in the ^/-direction but periodic 
parallel to the x-axis, we have 

cxt (r; w) = £ / — ^ e^^yy^Qn, q y , z; w) , (1) 

where q x , n = q x + nG, q n = yjql >n + q%, and 4> ext {q n , q y , z; u) satisfies ql4> ext = d 2 <fi cxt /dz 2 . 
Consequently, the general solution just outside (z ~ 0) or just within (~) the bi-layer system 
is given by 

where the upper (lower) sign in the exponential factor e ±QnZ may be chosen when z < or 
z > for the solution of Poisson's equation. The pair of coefficients T±(q n ;u) should be 
chosen to have different values above and below the plane at z = and are determined by the 



T + {q n -u)e qnZ + T_{q n -uj)e 
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boundary conditions, on the electrostatic potential and electric field Also, qii = (q x , q y ) is an 
in-plane 2D wave vector. The induced potential is obtained by solving V 2 0ind = — Pmd/eo- 
This may be rewritten as 



ind (r';u;) = / dv' v(v, v')p ind (v'; ui) 

= J dv' dv" u(r,r')x (0) (r', v";uj) <f) cxt (v" ; uj) , (2) 

where v(v, r') is the Coulomb interaction potential between the electrons in the 2D system, 
and x^( r '> r";w) is the nonlocal, frequency-dependent density-density response function. 
By Fourier transforming the response function and the external potential in the variables x 
and y, we obtain after a straightforward calculation the result 

ind (r; w) = - £ <? n (q||; w) e^ 2 e *.-*^w , (3) 

where the Fourier component of the surface response function is defined by 



dz' f dz" e q " z ' x (0) 0'> 2"; Qx,n, q y \ w)0 ext (g n , g„, 2") . (4) 



0« (q[|5 = - 

Combining the results in Eqs. (??) and (|3]), we obtain the total potential, which is the sum 
of the external and induced potential outside the layered system 

0(r; w) = J2 / 7^ [ e ~^ - 0n(q||5 w) e"" z ] e ^^ +i ^ for * £ , (5) 

n=-oo l^ 71 ") 

where q n , q X)Tl are defined above. We must now calculate the Fourier components of the 
surface response function <7 n (qii; uj) for a layered 2DEG. There is a dielectric medium with 
dielectric constant e filling the space between adjacent layers, except for a very thin vacuum 
region right next to the 2D layer. The electrostatic potential in the vacuum region at z = ze 
between layers at z = ze and z = z^+i is given by <pe(v, t) = e~ luJt <pe(v, ui), where 

J<(r; co) = J 7^1 [Ape-^*-** + B^W Z_ *<>] e iq ^ x+iq « y . (6) 



n=— oo 



Also, at the last interface we take 
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?L+i(r; uj) 



E 



rf 2 q|| 



ti+i e 



-q n {z-z L ) 



for z > z L+1 



(7) 



where the coefficient Vl+i must be determined from the matching boundary conditions like 
and However, we only choose the exponentially decaying term to ensure that the 
solution for the potential does not blow up. 

Both the potential <fi and its derivative edcf)/dz must be continuous at the vacuum- 
dielectric medium interface. However, the electric field is discontinuous across the 2D 
charged layer on which the induced electron sheet density is 



(n) 



(q n , u) = e\f\q n , u) (A^ + B^) , 



(8) 



where Xt fam w ) is a Fourier component of the single-particle density-density response func- 
tion of the E-th layer. We note that we employ no other parameter to describe the gated 
grating period. As a consequence, we only have to specify the grating potential in our 
numerical calculations by the period d. After some straightforward algebra, we obtain 



An) 
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(n) 



(n) 



Mi n) 



-g(q n , w) 



(9) 



where g(q n , oj) = g n (ci\\; u>) depends only on the total wave number q n = sjql~ n + q%, -M^ 
TW(ai) ® • • -® T(n)(ai) ® T^(a ), and 
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In Eq. (TTOl) . ae(q n , u) = (e 2 /2eoe) x\ (<7n> °°) is the polarization function for each 2D charged 
layer. When we equate the electrostatic potential just inside the material at the last layer 
to the electrostatic potential outside, we obtain from Eq. ffTUl) . after solving for the surface 
response function 
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[1 - e - 2ea L+l (q n , u)]M$ - [1 + e + 2ea L+1 {q n , u)\M<$ 



[l-e- 2ea L+1 {q n , u)]M® - [1 + e + 2ea L+1 (g n , w)]^? 

where we obtain the elements of the (2 x 2)-matrix M.^ 71 ' by evaluating the product of 
L transfer matrices T^ n ' whose elements are defined in Eq. fflOl) for a structure containing 
L + 1 layers. Clearly, in this formalism, the 2D charged layers are coupled through the 
Coulomb interaction and we assume that there is no interlayer electron hopping between 
layers. The surface response function in Eq. (TTTj) is a useful tool for calculating the normal 
mode spectrum of plasmon excitations for a finite number of layers and can also be employed 
to investigate the role played by layer separation on the loss function Im [g(q n , w)]. 14 The 
plasmon dispersion is obtained by setting the denominator in Eq. (11 ip equal to zero. That 
is, the plasmon resonances can occur for all values of n and are determined by the angle the 
in-plane polarization of the incident electromagnetic field makes with the x-axis. 

If there is a an electric field applied on the £-th layer creating a current and a drift velocity 
vd , then there is a Doppler shift on the angular frequency of the density-density response 
function of this layer and we have to replace u by u — q„ • v D , where q n = (q x , n , %)■ 

In the special case when there is a sing le layer, we set = M$ = 1 and M$ = 

M$ = in Eq. ([II]) which becomes^^ 

2 

Ssingleten, U)) = 1 - — — r . (12) 

1 + e + 2ea(q n , u) 

The surface response function for a semi-infinite slab of dielectric medium can then be 
deduced from this result by setting the polarization function a(q n , u) equal to zero. 



III. NUMERICAL RESULTS AND DISCUSSION 



In this section, we present the results of our numerical calculations of the real and imag- 
inary parts of the plasmon frequency for a bilayer 2DEG. Since inter-layer hopping is not 
included in our model, the layer separation must be chosen sufficiently large to satisfy this 
condition. This requires that the inter-layer separation is much larger than the lattice con- 
stant of the host material, which we assume to be GaAs/AlGaAs. In our calculations, we 
chose the separation between the two parallel layers of 2DEG to be a = 100A. The electron 
effective mass is chosen as 0.067 m e , appropriate for GaAs, where m e is the free-electron 
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mass, the electron density of each layer is taken as U2D = 2 x 10 11 cm~ 2 which is typical 
for a 2DEG. The corresponding Fermi energy and Fermi frequency are Ep = 7.14 meV and 
u>f = 10.85 s" 1 , respectively. We assume that the electrons in one layer have a drift velocity 
Vd = 2.5 Vp. The region between the two layers is filled with a dielectric material. The 
dielectric function of this material is approximated by the local frequency- dependent form, 

1 /2 

i.e., e{io) = 1 — (cup/o;) 2 , where u p = {AiinzDe 2 / e s ebm*) is the 3D plasma frequency of the 
dielectric material with background dielectric constant e& and e s = 4ireo€b. In our numerical 
calculations, the value that we used was u p = 0.92 Wi?. The whole system is embedded in a 
medium with dielectric constant 6^=13. 1. 

We justified using values for the drift velocity greater than the electron Fermi velocity by 
a factor of 2.5 by noting that we are in a region of Cherenkov-like instability of plasma waves 
where, the drift velocity may be greater than the phase velocity. Moreover, we have employed 
the force-balance equation^ as an alternative to the Boltzmann equation to calculate the 
drift velocity as a function of applied electric field and temperature in studying on nonlinear 
electron transport in impurity-limited quantum wires. This system is less computationally 
intensive but, nevertheless, gives a reasonable estimate of the drift velocity as a function of 
strong electric field and as well as at various temperatures. For example, when the ratio 
of the impurity to electron concentration is 0.01 and the temperature T m 100 K, the drift 
velocity can vary as 1.0 x 10 7 cm/s < v a < 3.0 x 10 7 cm/s for an electric field which varies 
over the range 0.5 V/m < E < 2.4 V/m. For the 2DEG density and the electron effective 
mass chosen, the value of the Fermi velocity is vf = 2.09 x 10 7 cm/s. Therefore, vd which 
we chose is reasonable in the frame of the force-balance formalism. 

We first examined the effect which we obtain on the complex plasmon solutions by varying 
ri3D- In Fig. 2, we plot the real (black dots) and imaginary (red dots) parts of the complex 
solutions u> (in units of up) of the plasmon dispersion equation in terms of the in-plane 
wave vector q» (in units of hp) for the parameters given above and for u p = 0.92^. No 
grating was applied in these calculations. Only those solutions with a nonzero imaginary 
part are shown in these plots. Then, in Fig. 3, we increased u p to twice its value, i.e., 
we chose u p = 1.84 ujp which corresponds to an increase of n 3 D by a factor of 4. All the 
other parameters were kept the same as in Fig. 2. We can see, by comparing Fig. 2 with 
Fig. 3, that 5m u in Fig. 3 has its local maximum at a higher value and exists over a wider 
range of wave vectors q». Specifically, in Fig. 2, the imaginary part of u has a maximum 
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at qu ~ 0.38 kp for which 5m u ~ 0.37^. In Fig. 3, we see that the maximum is located 
at 5|i ~ 0.62 kp where 3m u ~ 0.57a;i?. The increase in the value of u is by more 
than 50%. Therefore, the electronic collective excitations become more unstable, since their 
lifetime is inversely proportional to the imaginary part of the frequency. This enhancement 
of Qm u is reasonable since, in this case, an increase in n^n results in stronger interactions 
and stronger coupling between the two layers. We found three plasmon excitation modes, 
but in the figures we present two of them which have nonzero imaginary frequency. Both of 
them are acoustic, while the third mode which is not shown is optical. The mode which is 
the most unstable in both figures is the least energetic acoustic one. The mode just above 
this one, has a non-zero imaginary part too which is much smaller than the lowest one. The 
lowest plasmon mode bifurcates at q« fa 0.7kp. However, only the higher branch in this 
bifurcation has nonzero imaginary part and only this is shown. Compared with the model 
in Ref. 18), we include a spatial separation between the 2DEG layers and passed a current 



through one of them. This physical separation leads to the existence of the appearance of 
two unstable modes. The optical mode arises from the surface plasmon mode and is robust 
against excitations. This is unlike the acoustic modes which come about from the in-plane 
charge density fluctuations for which the screening properties are much different than the 
bulk. 

We next turn to the results of our investigation for a bilayer with a spatially modulated 
electron density. As described above in Sec. HI1 the modulation is applied along the in- 
direction. The period was taken to be d = 5 /im which corresponds to a reciprocal lattice 
vector G = 2ii/d « 1.26 x 10 6 m _1 . This chosen value for d is a typical experimental 
characteristic for metal gratings. The other parameters were chosen the same as in Fig. 2. 

We solved the dispersion equation g(q n ,u) = oo, where g is the surface response function 

r i V 2 

of the bilayer, q n = (q x + nG) 2 + q 2 ) is an in-plane wave vector and n is an integer. 
For n = 0, we have the case of the unmodulated bilayer which we mentioned above in 
Figs. 2 and 3. We are interested in the effect which the modulation has, i.e., for the case 
n = ±1, ±2, .... In Fig. 4 (a), we plot 3m u as a function of q x for three different values 
of n = 0, 5 and 10 when q y = 0. We see that as n assumes larger positive values, the 
corresponding 3m u curves are shifted to higher values of q x with negligible change in their 
local maximum value. In order to better understand the meaning of this shift, we examine 
a fixed q x = 0.6 kp- We then see that the curve that corresponds to n = 10 (blue dots) has 



8 



a much larger 3m to than the curves for n = (black dots) and n = 5 (red dots). Therefore, 
the presence of the modulation leads to a shifting which in turn produces an amplification of 
the instability for values of q x > 0.43 k F where the set of three curves corresponding to the 
acoustic branches intersect. However, for q x < 0.43 hp, we have the opposite effect leading 
to an attenuation. For negative values of n, we see from Fig. 4 (b), that the picture is more 
complicated since we have in addition to the g x -shift, a change in the shape of the curve 
for high values of n. As we can see though for low values of n, the result is the same as 
before, i.e., there are intervals of q x where we have amplification and intervals where we have 
attenuation. This is a result of the constructive or destructive mixing of the transmitted 
and reflected electromagnetic fields. This is determined by the wavelength of the applied 
field and, as a result, leads to some range of wave vectors when there is attenuation and 
another range where there is amplification. Comparing now Fig. 2 to Figs. 4 (a) and 4 (b), 
we can conclude that the presence of the modulation can give rise to an enhancement in the 
value of the 3m ui curves as a result of the stronger mixing of the incident and reflected 
electromagnetic waves. 

The effect on the imaginary part of the plasmon frequency u when we vary the period d 
of the modulation is shown in Figs. 5 (a) and 5 (b) for fixed n = 2 and n = —2, respectively. 
In Fig. 5 (a), we see that as the period of the modulation increases from d = 0.5 /xm (black 
dots) to d = 500 /iin (green dots), there are again intervals of q x where we have enhancement 
and intervals where we have attenuation. In addition, the curves approach a limiting value 
as the period d increases. The behavior is given approximately by the curve corresponding 
to d = 500 fim. This can be explained by calculating the reciprocal lattice vector G = 2n/d 
which is negligibly small as d becomes large. Large values of d correspond to the case that 
there is no modulation, or equivalently to the case that n = 0. Therefore, the limiting curve 
represents the response of the system when there is no spatial modulation from a grating. 
For n = —2, we see from Fig. 5 (b) that the picture is almost the same with the only 
difference being that the curve corresponding to d = 0.5 /xm is deformed from its parabolic 
shape. In both figures, there is an increase in the height of the curves when we increase 
the period from d = 0.5 fim to higher value. We emphasize that the significance of "n" is 
that it denotes the Bragg order numbers in the presence of the modulation. The optical 
polarization and the interference of a pair of optical-polarization waves are is, of course, 
determined by the value of n. This is why opposite values of n can lead to different results 
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since the polarization depends on the magnitude and sign of the Bragg number. 

IV. CONCLUDING REMARKS 

In this paper, we used a linear response function approach to investigate the plasmon 
instability for a bilayer system. In one of the 2DEG layers, a steady current is passed with 
drift velocity which induces a Doppler shift qv^ in the frequency-dependent polarization 
function. Here, we solved the plasmon dispersion relation in the complex frequency plane, 
following the work of Bakshi, et al.— ~— We extended their work by employing a model which 
couples the 2DEG layers to a bulk plasma which could be doped with free carriers and also 
by including the effect due to electron density modulation. The electrons in each 2D layer 
interact with the bulk plasma sandwiched between them. The crucial finding here is that 
there are two plasmon branches whose excitation frequency has nonzero imaginary part. 
They are both acoustic-like in nature. The lower-frequency mode is more unstable than the 
one with higher frequency. The optical mode, which is not shown in the figures, is a result 
of the hybridization between the surface plasmon and the 2DEG plasmon. We have also 
demonstrated how the instability of these modes is affected by varying the period of the 
modulation and the reciprocal lattice vector which couples the plasmon mode to an external 
electromagnetic field which is used to probe the system. 

The modulating grating on top of a conducting sheet generates and mixes Bragg modes 
of a reflected/transmitted electromagnetic field. The existence of the Bloch-like modes due 
to the modulation is a direct consequence of the nonlocal mixing of specular and diffrac- 
tion modes of the reflected electromagnetic field by free-electron induced optical polarization. 
Furthermore, there could be interference between a pair of surface optical-polarization waves 
with different Bragg order numbers in the presence of a modulation. The interference of 
these two counter-propagating surface waves leads to the formation of a Wannier-like state 
with associated electromagnetic fields localized within the gaps of the modulating poten- 
tial. These are the effects which contribute to the physical difference between the induced 
potential with and without a periodic potential. 
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2DEG 
2DEG 




FIG. 1: Schematic illustration of a bilayer with an applied electrostatic modulation periodic in the 
x direction. The surface at z = is covered by a grating (red lines) with period d and the width 
of the grating is assumed small compared to the period d. The doped barriers (AlGaAs) are the 
unshaded areas and quantum wells (GaAs, blue areas) in the z direction are indicated as the 2DEG 
layers. The system is embedded in a dielectric medium. 
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FIG. 2: Plots of the real (left scale) and imaginary (right scale) parts of the plasmon frequency for 
a bilayer 2DEG system with spacing a = 100 A between the layers. The frequencies are expressed 
in units of ojf and are plotted as functions of the in-plane wave vector tjii in units of kp. The 

real solutions, shown in black and labeled by ui^, and oo^\ have branches which bifurcate. The 
imaginary solutions, labeled correspondingly by u:^ and uj^ are shown in red. No modulating 
potential was applied in these calculations. Only those plasmon frequencies with nonzero imaginary 
parts are presented. 
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FIG. 4: (a) Imaginary part of the plasmon frequency in units of ujf for n = 10 (blue dots) has a 
much larger 9m ui than the curves for n = (black dots) and n = 5 (red dots). The parameters 
used in the calculations are the same as Fig. 2. The higher-frequency curves are for the acoustic- 
like plasmons discussed in Fig. 2. The group of curves with the lower frequency correspond to the 
optical plasmons. (b) The same as (a), except that n = 0, —5, —10. 
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FIG. 5: The imaginary part of the plasmon frequency u as a function of q x with q y = 0. We 
vary the period d of the modulation. Here, d = 0.5 //m (black dots), d = 50 /im (red dots) and 
d = 500 /im (green dots). In (a), n = 2 and (b) n = —2. All other parameters are the same as 
Fig. 2. 
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